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Measurement of the 20 and 90 keV resonances in the 18 0(p, a) 15 N reaction via THM 
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The ls O(p, a) 15 N reaction is of primary importance in several astrophysical scenarios, including 
fluorine nucleosynthesis inside AGB stars as well as oxygen and nitrogen isotopic ratios in meteorite 
grains. Thus the indirect measurement of the low energy region of the 18 0(p, a) 15 N reaction has been 
performed to reduce the nuclear uncertainty on theoretical predictions. In particular the strength 
of the 20 and 90 keV resonances has been deduced and the change in the reaction rate evaluated. 
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Understanding fluorine production would make predic- 
tions on asymptotic giant branch (AGB) star nucleosyn- 
thesis more accurate. The AGB stage represents the final 
nuclosynthetic phase in low and intermediate mass stars. 
AGB stars play an extremely important role in astro- 
physics because heavy elements along the stability valley 
are produced in their interiors through slow neutron cap- 
tures (s-process). 19 F can be produced during the ther- 
mal pulse that is ignited in the 4 Hc-rich intcrshcll region, 
following the ingestion of the 13 C pocket. The subsequent 
third dredge-up (TDU) episode mixes the products of 
shell flash He-burning, including fluorine, and s-process 
nuclei to the outer layers. Therefore 19 F abundance con- 
stitutes a key parameter to constrain AGB star models 
[l[ . Anyway, if the theoretical abundances are compared 
to the observed ones, a remarkable discrepancy shows up 
because the largest 19 F abundances cannot be matched 
for the typical 12 C/ 16 ratios [lj. A possible way to 
explain 19 F abundance and several related observables 
(such as isotopic ratios in meteorite grains Q) can be 
linked to 18 as the 18 0(p, a) 15 N reaction represents the 
main 15 N production channel, both in the intershell re- 
gion and at the bottom of the convective envelope [l], Q • 
During the thermal pulse 15 N is burnt to 19 F via the 
15 N(a,7) 19 F reaction. Thus a larger 18 0(p, a) 15 N reac- 
tion rate would lead to an increase of the 19 F supply as 
well as to an enrichement of 15 N in the stellar surface, 
which would play a key role to explain the long-standing 
problem of the 14 N/ 15 N ratio in meteorite grains (Q and 
references therein). 

The 18 0(p, a) 15 N reaction has been the subject of sev- 
eral experimental investigations 0, 0] and many features 
are known from spectroscopic studies [!, H, 0, H[ • Never- 
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theless the reaction rate for the process has a consider- 
able uncertainty [9]. In the — 1000 keV energy range, 
where 9 resonances occur, the reaction rate is essentially 
determined by the 20 keV, 144 keV and the 656 keV res- 
onances [9j. With regard to the 20 keV resonance, its 
strength is known only from spectroscopic measurements 
performed through the transfer reaction 18 0( 3 He, <i) 19 F 
Q and the direct capture reaction 18 0(p, 7) 19 F Such 
estimates, which are based on the deduced spectroscopic 
factors, are strongly model dependent (being connected 
to the adopted optical model potentials) and affected by 
large and not-well-defined uncertainties. An additional 
important source of uncertainty on the reaction rate is 
connected with the determination of the resonance en- 
ergy [g]. Furthermore the spin and parity of the 8.084 
MeV level in 19 F (corresponding to a 90 keV resonance 
in the 18 0(p, a) 15 N cross section) has not been estab- 
lished. The uncertainties on nuclear physics inputs have 
made astrophysical predictions far from conclusive [|J ■ 

In order to reduce the nuclear uncertainties affecting 
the reaction rate estimate we have performed an exper- 
imental study of the 18 0(p, a) 15 N reaction by means of 
the Trojan horse method (THM), which is an indirect 
technique to measure the relative energy-dependence of a 
charged-particle reaction cross section at energies well be- 
low the Coulomb barrier ([I3,[ll[ and references therein). 
The cross section of the relevant A + x — > c + C reaction 
is deduced from a suitable A + a(x © s) — » c + C + s pro- 
cess, performed in quasi-free (QF) kinematics. The beam 
energy is chosen larger than the Coulomb barrier for the 
interacting nuclei, so the break-up of nucleus a (the so- 
called Trojan-horse) takes place inside the A nuclear 
field. Therefore, the cross section of the A + x — > c + C 
reaction is not suppressed by the Coulomb interaction of 
the target-projectile system. 

In a previous investigation, carried out at the Cy- 
clotron Institute, Texas A&M University, Texas (USA) 
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[T3 |. the 18 0(p, a) 15 N reaction has been measured via 
the THM through the 2 H( 18 0, a 15 K)n process in the 
— 1000 keV 18 — p relative energy range. For the 
first time the energy region below 70 keV had been in- 
vestigated. In the present work we focus on a new study 
of the 18 0(p, ct) 15 N reaction by means of the same THM 
process. The aim is to span the 18 0— p relative energy re- 
gion below 250 keV with an improved energy resolution, 
in order to deduce resonance parameters and 3 n values 
of the 8.014 and 8.084 MeV 19 F levels. The experiment 
was performed at Laboratori Nazionali del Sud, Cata- 
nia (Italy). The SMP Tandem Van de Graaff accelerator 
provided the 54 MeV 18 O beam impinging onto thin self- 
supported deuterated polyethylene (CD2) targets. The 
detection setup consisted of a telescope (A), devoted to 
15 N detection, made up of an ionization chamber and a 
silicon position sensitive detector (PSD A) on one side 
with respect to the beam direction, and three additional 
silicon PSD's (B, C and D) on the opposite side. An- 
gular conditions were selected in order to maximize the 
expected QF contribution. Channel selection has been 
accomplished by gating on the AE — E 2D spectra to 
select the nitrogen locus. 

Compelling evidence for the occurrence of the 
QF mechanism is given by the shape of measured 
momentum-distribution, if it follows the shape of the 
deuteron wave function. In the analysis, the theoreti- 
cal distribution, given by the square of the Hulthen wave 
function in momentum space in the plane wave approxi- 
mation [111], is superimposed onto the experimental one. 
The good agreement demonstrates that the QF mecha- 
nism is present and dominant in the \p$\ < 50 MeV/c 
neutron momentum range. Thus, in the following anal- 
ysis only the phase space region for which the \p 3 \ < 50 
MeV/c condition is satisfied is taken into account. 

Angular distributions for the 18 0(p, a) 15 N reaction 
were extracted as discussed in [ll|. They were used to 
perform the necessary validity test on the deduced cross 
section and to infer the spin-parity for the states. The 
a emission angles in the cm. system (9 c . m .), covered 
in the present experiment, were about Q c .m. = 0° - 60°, 
0c. m . = 40° - 110° and 6 c . m . = 90° - 150° for the A- 
B, A-C and the A-D detector coincidences respectively. 
The presence of an overlap region allowed for relative 
normalization between the cross sections deduced from 
each couple. Angular distributions were extracted for 
several energies, focusing in particular on the 20, 90 and 
143.5 keV resonances. The results are displayed in Fig. 
Q] where the experimental data are given as filled circles 
(20 keV), squares (90 keV) and triangles (143.5 keV). 
Errors on the half-off-energy- shell (HOES) [llj cross sec- 
tion 4S account for statistics and for the deconvolution 
of the single resonance contributions. The error on c , m , 
represents the width of each bin (chosen in order to have 
reasonable statistical precision). 

From Fig. Q]it turns out that the = \ + assignment 
for the 143.5 keV resonance is confirmed, the angular dis- 
tribution for that level being isotropic (see Fig. [T|). This 
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FIG. 1: Experimental angular distributions for the 
18 0(p, a) 15 N reaction for the three resonances in the — 250 
keV energy range. The full lines come from the fitting with 
the curve of Eq. [1] 



result represents a cross check of the method, since we 
are able to reproduce the angular distribution for a well 
known resonance. Therefore we extracted the angular 
distribution for both the 90 keV and 20 keV resonances. 
A fit of the experimental data was performed by: 
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where L is the smallest among the spin J of the formed 
compound nucleus and the angular momentum I quan- 
tum numbers in the entrance and the exit channels [l3|. 
The best fit for the 90 keV resonance is achieved for L = 1 
(x 2 = 0.67), thus we infer a spin for the 8.084 MeV ex- 
cited state of 19 F of |. We consider that this L value 
corresponds to the angular momentum in the exit chan- 
nel, because of considerations on the proton width of this 
level, and the parity to be positive (the 15 N ground state 
being | ). For the 20 keV resonance, the best fit is ob- 
tained for L — 2 (x 2 = 3.11), supporting the J w = | + 
spin-parity assignment for the 19 F excited state at 8.014 
MeV (see @, Q). In the extraction of the angular dis- 
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FIG. 2: Cross section of the TH reaction (full circles). (1), 
(2) and (3) refer to the 8.014, 8.084 and 8.138 MeV 19 F lev- 
els. The full line represent the result of a fit including three 
Gaussian curves and a 1st order polynomial. 



tributions the 18 O — p relative energy is kept fixed thus 
preventing possible bias caused by the HOES nature of 
the deduced cross section. The HOES effect might be ex- 
pressed by a constant renormalization factor, and there- 
fore the angular distributions are given in arbitrary units 
in Fig. CD 

The extracted HOES differential cross section has been 
integrated in the whole angular range. It was assumed 
that in the region where no experimental angular distri- 
butions are available, their trend is given by the fit of Eq. 
[TJ The resulting 2 H( 18 0, a 15 N)n reaction cross section is 
shown in Fig. [2] (full circles). The experimental energy 
resolution turned out to be about 40 keV (FWHM), in 
agreement with the value predicted by Monte Carlo sim- 
ulations. Horizontal error bars represent the integration 
bin while the vertical ones arise from statistical uncer- 
tainty and angular distribution integration. Indeed, be- 
cause of the limited energy resolution, resonances were 
not well resolved thus each contribution to the reaction 
yield had to be disentangled when extracting angular dis- 
tributions. The solid line in the figure is the sum of three 
Gaussian functions to fit the resonant behavior and a 
straight line to account for the non-resonant contribution 
to the cross section. This fit was performed with the sole 
aim of extracting the resonance energies: E^ = 19.5±1.1 
keV, E fla = 96.6 ± 2.2 keV and E Rs = 145.5 ± 0.6 keV 
(in fair agreement with the ones reported in the litera- 
ture [§]) and to deduce the peak values of each resonance: 
iVi = 138 ± 8, N 2 = 82 ± 9 and N 3 = 347 ± 8 (arbitrary 
units) . The peak values were used to derive the resonance 
strengths that are the relevant parameters for astro- 
physical application in the case of narrow resonances Q . 
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FIG. 3: Ratio of the THM reaction rate (full black line) to the 
NACRE one [0]. For comparison, the NACRE rate is shown 
as a red line. The dot-dashed and dotted lines show the upper 
and lower limits allowed by experimental uncertainties. 



The THM cross section for the A + a(x © s) — > c + C + s 
reaction proceeding through a resonance Fi in the sub- 
system F = A + x = C + c is [U G3] 



dE Cc dfl s (E - E Ri y + Tf{E)/A- 



(2) 



Here, Mi(E) is the direct transfer reaction amplitude for 
the binary reaction A + a — > F^ + s populating the res- 
onant state Fi with the resonance energy E Ri , E is the 
A — x relative kinetic energy related to Eq c by the en- 
ergy conservation law, T^ Cc ^.{E) is the partial resonance 
width for the decay Fi — ► C + c and Tj is the total reso- 
nance width of Fi . The appearance of the transfer reac- 
tion amplitude Mi(E) instead of the entry channel partial 
resonance width T^ x ^ . (E) is the main difference between 
the THM cross section and the cross section for the reso- 
nant binary sub-reaction A + x — > C + c [ll], Q3] • On the 
other hand, the resonance parameters deduced by means 
of the THM are not affected by the electron screening, 
distorting the direct S-factor at the astrophysically rele- 
vant energies [lj|. 

The peak THM cross section taken at the E^ reso- 
nance energy for the (p, a) reaction A + x — >C + c is 
given by 



Ni = A 



T ai (E Rt )M?(E Ri ) 
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where Tc Ci (E) = T ai (E) . In this work we did not ex- 
tract the absolute value of the cross section. The proton 
and alpha partial widths for the third resonance at 143.5 
keV are well known Q , thus we can determine the alpha 
partial width and the strength for the first and second 
resonances from the ratio of the peak values of the THM 
cross sections: 
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Here, (w 7 )< = w; T^S^) Tp^^/r^J is the i-th 
resonance strength, Wi = J%/(Ja Jx) is the statistical fac- 
tor, J = 2 J + 1, Ji is the spin of the i-th resonance 
and Ja, J x are the spins of A and x, respectively. When 
determining (a> 7 )j the effect of energy resolution in our 
experiment has been taken into account. The electron 
screening gives a negligible contribution around 144 keV 
(4% maximum (l5j). thus no systematic uncertainty is 
introduced by normalizing to the third resonance. In 
the plane wave approximation Mj ks ip a (p sx ) Wa x (pax), 
where <p a (Psx) is the Fourier transform of the s-wave ra- 
dial s—x bound wave function, p sx is the s—x relative mo- 
mentum, and Wax = < Iax \Vax\pAx > is the form factor 
for the synthesis A + x — > F{ , I^ x is the overlap function 
of the bound state wave functions of A, x and the reso- 
nant wave function of F* , pa x is the A — x relative mo- 
mentum. In practical calculation we approximated /jT*, 

1 /2 

by S i ( P(Ax) i j where Si is the spectroscopic factor for 
configuration Fi — A + x and ip(Ax)i is the single-particle 
bound-state type wave function describing the resonance 

r-n r (E ) 

state Fi. Since in Eq. 2] only the ng^jg^ra ratios show 
up, the dependence on spectroscopic factors and pro- 
ton widths is completely removed and, as a consequence, 
(u)j)i is connected to (077)3 through the easily calcula- 
ble ratios of the single-particlejwidths to the form factors 
Wax- If (077)3 is taken from [161 ] ■ by means of Eq. [4] one 
gets {oj^) l = 8.3tl 8 6 x 10 ~ 19 eV ' which is well within the 
upper and lower limits given by NACRE, 6+* 7 x 10" 19 eV 
[9j . While NACRE recommended value is based on var- 
ius kinds of estimates, the present result is obtained from 
experimental data, thus the accuracy of the deduced res- 
onance strength has been greatly enhanced. The largest 
contribution to the error is due to the uncertainty on 
the resonance energy, while statistical and normalization 
errors sum up to about 9.5%. With a similar approach 
we have obtained (077)2 = (1.76 ± 0.33) x 10~ 7 eV (sta- 
tistical and normalization errors ~ 13%) for the 90 keV 
resonance, in good agreement with the strength given by 
NACRE, (1.6 ± 0.5) x 10~ 7 eV @. This gives us confi- 
dence in the theory used in the THM allowing for a cross 
check of the method. 

By using the narrow resonance approximation Q, 
which is fulfilled for the resonances under investigation, 
the reaction rate has been deduced and compared with 



the one reported in NACRE @. In Fig. [3] the ratio 
of the THM reaction rate to the NACRE one for the 
ls O(p, a) 15 N is shown as a full black line. The dot-dashed 
and dotted black lines represent the upper and lower lim- 
its respectively, allowed by the experimental uncertain- 
ties. In the low temperature region (below T9 = 0.03, 
Fig. [3ji) the reaction rate can be about 35% larger than 
the one given by NACRE (full red line) , while the indeter- 
mination is greatly reduced with respect to the NACRE 
one (dot-dashed and dotted red lines mark the upper and 
lower limits, Fig. [3]). Those temperatures are typical of 
the bottom of the convective envelope, thus an increase 
of this reaction rate might have important consequences 
on the cool bottom process [2] and, in turn, on the sur- 
face abundances and isotopic ratios in AGB stars. The 
8.084 MeV excited state of 19 F (corresponding to the 90 
keV resonance) provides a negligible contribution to the 
reaction rate in agreement with previous estimate 0]. 

In conclusion, in this paper we demonstrated, for the 
first time, the power of indirect THM which allowed to 
determine the strength of the low-lying 20-keV resonance 
in 19 F, elusive for any direct technique. This resonance 
plays an important role in the determination of the reac- 
tion rates of the key astrophysical reaction 18 0(p, a) 15 N. 
We also presented a new way of determination of the res- 
onance strength by measuring the relative strengths of 
the known and unknown resonances and using the half- 
off-sell R-matrix determine the strength of the unknown 
resonance avoiding information about the spectroscopic 
factors. The parameters of the 20 keV resonance in the 
18 0(p, cv) 15 N reaction relevant for astrophysics together 
with those for the nearby resonances could be extracted 
and importantly their uncertainty was strongly reduced 
(by a factor ~ 8.5), thanks to the newly developed ap- 
proach, which is based on experimental data in contrast 
to the NACRE one that relies on varius kinds of estimate. 
Indeed, in Eq. 0] only the ratio of the model dependent 
parameters in Eq. Oshows up, thus systematic uncertain- 
ties cancel out. In addition, our results are not affected 
by the electron screening, which can enhance the cross 
section by a factor > 2.4 at 20 keV [15| . thus spoiling 
any direct measurement of this resonance. At higher tem- 
peratures, higher energy resonances in the ls O(p, a) 15 N 
reaction can play a role. Their study will be the subject 
of forthcoming works. 
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